MLABS.021A PATENT 
SYSTEMS AND METHODS FOR DETERMINING BLOOD OXYGEN 
SATURATION VALUES USING COMPLEX NUMBER ENCODING 



REFERENCE TO RELATED APPLICATIONS 
[0001] The present application claims priority benefit under 35 U.S.C. §1 19(e) 
from U.S. Provisional Application No. 60/430,834, filed December 4, 2002, entitled, 
"Systems and Methods for Determining Blood Oxygen Saturation Values Using Complex 
Number Encoding" which is incorporated herein by reference. 

FIELD OF THE INVENTION 
[0002] The present invention relates to the field of pulse photometery. More 
specifically, the invention relates to calculating continuous saturation values using complex 
number analysis. 

BACKGROUND OF THE INVENTION 
[0003] Pulse photometry is a noninvasive technique for measuring blood analytes 
in living tissue. In this technique, multiple light sources emit light of differing wavelengths, 
which is transmitted through or reflected from a vascular bed. One or more photodetectors 
then detect the transmitted or reflected light as an optical signal. As the photons propagate 
through the tissue, they are subjected to random absorption and scattering processes due to 
the nonhomogeneous nature of the tissue. These effects manifest themselves as a loss of 
energy in the optical signal, and are generally referred to as bulk loss. In addition to bulk 
loss, the optical signal is modulated by the flow of arterial blood into the vascular bed. 
Moreover, the movement of venous blood into or out of the tissue, local tissue compression 
and local muscle movements super-impose yet another modulation on the optical signal, 
usually of lower frequency than the arterial flow. For example, Figure 1 illustrates detected 
optical signals that include the foregoing attenuation, arterial flow modulation, and low 
frequency modulation. Each optical signal, with its combined attenuation and modulations, 
e.g., each combined optical signal, is generally referred to as a photoplethysmograph 
(photopleth.) 
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[0004] Pulse oximetry is a special case of pulse photometry where the 
oxygenation of arterial blood is sought in order to estimate the state of oxygen exchange in 
the body. In order to calculate the oxygen saturation of arterial blood, two wavelengths of 
light, e.g. Red, at about 660 nm, and Infrared, at about 900 nm, are used to calculate the ratio 
of two dominant hemoglobin components, oxygenated hemoglobin (HB0 2 ) and 
deoxygenated hemoglobin (HB). The detected optical signals, which correspond to the Red 
and Infrared wavelengths, are first normalized in order to balance the effects of unknown 
source intensity as well as unknown bulk loss at each wavelength. The arterial pulses are 
then isolated by filtering each normalized signal, where a high pass or a band-pass filter takes 
advantage of the typically higher frequency of the pulsatile arterial blood, hence the name 
pulse oximetry. This normalized and filtered signal is referred to as the AC component and 
is typically sampled with the help of an analog to digital converter with a rate of about 30 to 
about 100 samples/second. For example, Figure 2 illustrates the optical signals of Figure 1 
after they have been normalized and bandpassed. 

[0005] In order to estimate blood oxygenation, a (Red/Infrared) ratio is calculated 
by dividing the strength of the Red AC (RdAC) by the corresponding strength of the Infrared 
AC (IrAC). The (RdAC/IrAC) ratio is then generally plugged into an empirical calibration 
curve equation that relates it to blood oxygenation. For example, reference can be made to 
Japanese Patent No. Sho 50/1975-128387, issued to Aoyagi, entitled "Optical Type Blood 
Measuring Equipment" 

[0006] The arterial blood flow generally has a higher fundamental frequency than 
other components of the photopleth, however, there are cases where the two frequencies may 
overlap. One such example is the effect of motion artifacts on the optical signal, which is 
described in detail in U.S. Patent No 6,157,850, issued to Diab et al., entitled "Signal 
Processing Apparatus" Another effect occurs whenever the venous component of the blood 
is strongly coupled, mechanically, with the arterial component. This condition leads to a 
venous modulation of the optical signal that has the same or similar frequency as the arterial 
one. Such conditions are generally difficult to effectively process because of the overlapping 
effects. 
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[0007] As described in the Aoyagi patent, the strength of each AC waveform may 
be estimated by measuring its size through, for example, a peak-to-valley subtraction, by a 
root mean square (RMS) calculations, integrating the area under the waveform, or the like. 
These calculations are generally least averaged over one or more arterial pulses. It is 
desirable, however, to calculate instantaneous ratios (RdAC/IrAC) that can be mapped into 
corresponding instantaneous saturation values, based on the sampling rate of the photopleth. 
However, such calculations are problematic as the AC signal nears a zero-crossing where the 
signal to noise ratio (SNR) drops significantly. For example, dividing two signals with low 
SNR values can render the calculated ratio unreliable, or worse, can render the calculated 
ratio undefined, such as when a near zero-crossing area causes division by or near zero. To 
try to avoid division by zero, the Ohmeda Biox pulse oximeter calculated the small changes 
between consecutive sampling points of each photopleth in order to get instantaneous 
saturation values. Figure 3 illustrates various techniques used to try to avoid the foregoing 
drawbacks related to zero or near zero-crossing, including the differential technique 
attempted by the Ohmeda Biox. 

[0008] Note that Ohmeda' s differential technique is equivalent to a calculation 
over a derivative of the photopleth, and the derivative has the same low SNR problem 
whenever a flattened section of the photopleth is used in the ratios calculations. For example, 
the derivative will have a zero or near zero value and the RdAC/IrAC ratio will become 
unreliable or undefined, even in a substantially noise free signal. For example, Figure 4 
illustrates the derivative of the IrAC photopleth plotted along with the photopleth itself. As 
shown in Figure 4, the derivative is even more prone to zero-crossing than the original 
photopleth as it crosses the zero line more often. Also, as mentioned, the derivative of a 
signal is often very sensitive to electronic noise. For example, according to "Pulse Oximetry: 
Analysis of Theory, Technology, and Practice" Journal of Clinical Monitoring, Vol. 4, Oct 
1988, a published paper by the designers of the Ohmeda Biox, the calculated instantaneous 
saturations over some sections of the photopleth can be off by more than 50 percent (0/0) 
from the real value over a time as short as 1/1 0th of a second. As the designers described in 
their paper, this result is clearly an artifact of the signal processing technique employed in the 
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Biox pulse oximeter since the blood saturation value can not change by that amount in 1/1 Oth 
of a second. 

[0009] Because of some of the foregoing drawbacks associated with the 
determination of instantaneous or point-by-point saturation from RdAC/IrAC ratios, 
designers now typically unequally weigh the calculated instantaneous saturation values over 
each photopleth, even when the photopleth is substantially noise free, with the consequence 
that a significant number of saturation values receive insignificant weights. This is 
tantamount to filtering out or ignoring valid signal data during the troublesome sections 
described above. 

SUMMARY OF THE INVENTION 

[0010] The result of the foregoing drawbacks is that the previous attempts fail to 
determine a stream of point-by-point saturation values. However, it is noteworthy that the 
sensitivity of the ratio calculation near a zero-crossing section of the waveform is not 
intrinsic to the photopleth itself, but rather an artifact of passing the detected signal through a 
high-pass filter. Accordingly, an aspect of the present invention includes a method of 
determining continuous and reliable calculations of the (RdAC/IrAC) ratio for each sampling 
point without concern for zero-crossing areas. As discussed in the foregoing and disclosed in 
the following, such determination of continuous ratios is very advantageous, especially in 
cases of venous pulsation, intermittent motion artifacts, and the like. Moreover, such 
determination is advantageous for its sheer diagnostic value. 

[0011] For purposes of summarizing the invention, certain aspects, advantages 
and novel features of the invention have been described herein. Of course, it is to be 
understood that not necessarily all such aspects, advantages or features will be embodied in 
any particular embodiment of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0012] A general architecture that implements the various features of the invention 
will now be described with reference to the drawings. The drawings and the associated 
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descriptions are provided to illustrate embodiments of the invention and not to limit the scope 
of the invention. 

[0013] Figure 1 illustrates a photopleths including detected Red and Infrared 

signals. 

[0014] Figure 2 illustrates the photopleths of Figure 1, after it has been 
normalized and bandpassed. 

[0015] Figure 3 illustrates conventional techniques for calculating strength of one 
of the photopleths of Figure 2. 

[0016] Figure 4 illustrates the IrAC photopleth of Figure 2 and its derivative. 

[0017] Figure 4A illustrates the photopleth of Figure 1 and its Hilbert transform, 
according to an embodiment of the invention. 

[0018] Figure 5 illustrates a block diagram of a complex photopleth generator, 
according to an embodiment of the invention. 

[0019] Figure 5A illustrates a block diagram of a complex maker of the generator 
of Figure 5. 

[0020] Figure 6 illustrates a polar plot of the complex photopleths of Figure 5. 
[0021] Figure 7 illustrates an area calculation of the complex photopleths of 
Figure 5. 

[0022] Figure 8 illustrates a block diagram of another complex photopleth 
generator, according to another embodiment of the invention. 

[0023] Figure 9 illustrates a polar plot of the complex photopleth of Figure 8. 

[0024] Figure 10 illustrates a three-dimensional polar plot of the complex 
photopleth of Figure 8. 

[0025] Figure 11 illustrates a block diagram of a complex ratio generator, 
according to another embodiment of the invention. 

[0026] Figure 12 illustrates complex ratios for the type A complex signals 
illustrated in Figure 6. 

[0027] Figure 13 illustrates complex ratios for the type B complex signals 
illustrated in Figure 9. 
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[0028] Figure 14 illustrates the complex ratios of Figure 13 in three (3) 
dimensions. 

[0029] Figure 15 illustrates a block diagram of a complex correlation generator, 
according to another embodiment of the invention. 

[0030] Figure 16 illustrates complex ratios generated by the complex ratio 
generator of Figure 1 1 using the complex signals generated by the generator of Figure 8. 

[0031]. Figure 17 illustrates complex correlations generated by the complex 
correlation generator of Figure 15. 

[0032] Figure 18 illustrates the square root of the magnitude of the complex ratios 
of Figure 16 vs. the complex correlations of Figure 17. 

[0033] Figure 19 illustrates a plot of the instantaneous saturation of the data used 
to generate the complex ratios shown in Figure 12, as well as the corresponding complex 
photopleth from which the saturation was calculated. 

[0034] Figure 20 illustrates an expanded view of the saturation distribution results 
of Figure 19. 

[0035] Figures 21 and 22 illustrate Infrared and Red photopleths, respectively, 
modulated by venous pulsation. 

[0036] Figure 23 illustrates an instantaneous saturation which highlights a large 
spread of values over one photopleth, as compared to that of Figure 19 and Figure 20. 

[0037] Figure 24 illustrates a histogram of the distribution of the instantaneous 
saturation values of Figure 23. 

[0038] Figure 25 and 26 illustrate photopleth signals that are corrupted by motion 
artifacts. 

[0039] Figure 27 illustrates a polar plot for type A complex waveforms generated 
using the photopleths of Figures 25 and 26. 

[0040] Figure 28 illustrates a polar plot for type B complex waveforms generated 
using the photopleths of Figures 25 and 26. 

[0041] Figure 29 illustrates the complex ratios of Figure 27 after being filtered, 
according to embodiments of the invention. 
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[0042] Figure 30 illustrates a histogram of the filtered saturation points of 
Figure 29. 

[0043] Figure 31 illustrates plots of unfiltered instantaneous saturation values 
generated from the photopleths of Figures 25 and 26, as well as the phase filtered saturation 
values of Figures 29 and 30. 

[0044] Figure 32 illustrates magnitudes of complex frequency ratios calculated 
from the fundamental and harmonics of the photopleths of Figure 2. 

[0045] Figure 33 illustrates frequency transformed photopleths and phases of the 
corresponding ratios. 

[0046] Figure 34 illustrates the polar plot of the complex frequency ratios of 
Figure 32 and Figure 33. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0047] Complex numbers were generally invented/discovered in the 1500 f s by 
Galiermo Cardano in Italy as he was struggling to solve a general third order polynomial 
equation. Later, Argand suggested that each complex number may be represented as a point 
in a plane where its imaginary part is plotted on the Y-axis and its real part on the X-axis, this 
is referred to as the Cartesian form of a complex number. An alternate representation of 
complex numbers is called the polar form where a magnitude and an angle can designate a 
unique point in the Argand plane, hence representing a complex number. Although 
mathematicians at that time looked at complex numbers with great suspicion, it turned out 
that they were useful for the general solutions of polynomial equations and applicable in such 
diverse fields as quantum mechanics to describe the state of elementary particles. 

[0048] In the field of electrical engineering, students learn that circuit analysis, 
under alternating voltage conditions, can be significantly simplified if the concept of complex 
currents and voltages is introduced. When a voltage is applied to an electrical element, a 
current is caused to pass through it. Dividing the applied voltage by the corresponding 
current gives the resistance of the electrical element. In the complex domain, dividing the 
complex voltage by the complex current give rise ' to a complex form of resistance called 
"impedance," which is represented again by a complex number. The real part of the 
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impedance is the resistance while the complex part is related to capacitance and inductance. 
The complex part affects the phase, lead or lag, of the driving signals. It can be shown that 
the response of a linear system to a sinusoidal excitation is a sinusoid with the same 
frequency but generally of different amplitude and phase. A resistor element, for example, 
affects the amplitude only, whereas a capacitor or an inductor affects the phase. A 
combination of resistors, capacitors and/or inductors can affect the amplitude as well as the 
phase of the driving excitation and are considered examples of a linear system. Since each 
complex number consists of amplitude and phase, it is natural that they be used to encode the 
amplitude and phase and track their evolution throughout a linear system. Therefore, when a 
linear system is driven by a complex excitation, its output will be the same complex input 
multiplied by a complex scaling factor that scales its amplitude appropriately and adds a 
certain phase in accordance with the rules of complex multiplication. Note that the phase in 
this context is a relative one between two fundamental variables of the system, e.g. for an 
electrical circuit they might be voltage and current at a certain node of the circuit, or output 
voltage and input voltage, or the like. 

[0049] In pulse oximetry, there is no direct analog for voltages or currents. 
Generally, pulse oximetry deals with two highly correlated optical signals, e.g., the Red and 
Infrared signals, with fundamentally little or no discernible phase difference. Thus, encoding 
those signals as complex numbers seems to add little or no value to the signal processing. 
However, there are several conditions under which a variable phase difference may be 
introduced between the Red and Infrared signals. For example, motion artifacts create a 
condition where the sensor may decouple from the skin. In such a condition, the detected 
optical signals will have components that depend on the refraction through the sensor 
material itself instead of the wavelength of light, as well as the desired components that have 
traveled through the vascular bed. Venous pulsation creates another condition, which, as 
disclosed in the foregoing, affects the phase difference. Under the foregoing conditions 
where the phase is changing, use of complex number encoding provides advantages in signal 
processing, including providing the ability to continuously monitor the arterial saturation vs. 
time without concern about signal zero-crossing, as disclosed in the following. 
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[0050] A Hilbert Transformer is a signal processing technique that takes a real 
signal and converts it into a related signal which has its frequency components shifted by n/2 
radians for positive frequencies, and by -7i/2 radians for negative frequencies, without 
affecting their respective amplitudes. The book, "Theory and Application of Digital Signal 
Processing, Prentice-Hall f Inc." by Rabiner and Gold, introduces the subject. Formation of 
a complex signal is accomplished by considering the original signal itself as the real part of a 
complex signal and the output of the Hilbert Transformer as the imaginary part of the same 
complex signal. Such signals are generally referred to as an "analytical signal" in the signal 
processing field because the magnitudes of its negative frequencies are equal to zero. In the 
context of pulse photometry, we shall refer to such a complex signal as complex photopleth 
or complex AG. 

[0051]. It is noteworthy that the prior methodology of seeking a derivative of a 
signal shifts that signal's components by n/2 radians, similar to the Hilbert transform. 
However, the derivative of a signal also multiplies each corresponding amplitude by the value 
of the radian frequency co, thus magnifying the signal's frequency components as the 
components increase in frequency. Thus, the derivative is often more sensitive to electronic 
noise. In contrast, the Hilbert Transformer generally has a flat response with respect to 
frequency. However, a skilled artisan will recognize other transforms, derivatives, or the 
like, can be used to encode the imaginary part of the complex photopleth. 

[0052] When the foregoing complex encoding is applied to both RdAC and IrAC 
photopleths, as shown in a complex photopleth generator of Figure 5, then both complex 
RdAC and IrAC are in phase. As shown in Figure 5, the complex photopleth generator 500 
includes one or more log filters 502, one or more high pass filters 504, one or more Hilbert 
transformers 506, and one or more complex makers 508. The log and high pass filter 502 
and 508 generally normalize and filter the signals, as disclosed in the foregoing with 
reference to Figure 2. The Hilbert transformers 506 converts the real signal into related 
signals shifted by n/2 and -n/2 radians, also as disclosed in the foregoing. The complex 
makers 508 combine the output of the Hilbert transformers 506 with the input of the Hilbert 
transformers 506 to generate complex in phase photopleths, generally referred to herein as 
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type "A" complex signals, as shown in Figure 5A. Figures 6 and 7 illustrate polar plots of 
the complex photopleths generated from the generator 500 of Figure 5. 

[0053] Figure 8 illustrates a complex photopleth generator 800, according to 
another embodiment of the invention. As shown in Figure 8, the generator 800 includes the 
one or more log filters 502, the one or more high pass filters 504, the one or more Hilbert 
transformers 506, and the one or more complex makers 508 disclosed with respect to 
Figures. However, in the generator 800, the Hilbert transformer 506 accepting the 
exemplary Red signal encodes the real components of the RdAC signal rather than the 
complex components. Thus, the generator 800 of Figure 8 generates a complex photopleth, 
generally referred to herein as type "B" complex signals. Figures 9 and 10 illustrate polar 
plots of the complex photopleths generated from the generator 800 of Figure 8. Note that as 
time progresses, the locus of the complex points rotates around the origin but does not pass 
through it, regardless of type A or type B complex signals. Accordingly, the strength of each 
of the foregoing AC photopleth complex signals can be encoded in their respective 
magnitudes, which is the length between the origin and any point on the complex waveform 
in Figures 6, 7, 9 and 10. The morphology of the complex signals depend on the condition of 
the subject, such as, for example, age, blood pressure, arterial impedance, posture, or the like. 

[0054] Based on the foregoing disclosure, the complex photopleths of Figures 6, 
7, 9 and 10 can be classified with techniques like the "Slant Line Transform" or other 
classification techniques available in the field of image recognition or the like! Such 
classification techniques can be advantageously employed to help reject photopleths that are 
corrupted by noise, which can be an invaluable during episodes of motion artifacts. 

Generating Complex Ratios 

[0055] Figure 11 illustrates a complex ratio generator 1100 for generating what 
will generally be referred to as the "complex ratio." The complex ratio is a point-wise 
complex division of type A or type B RdAC and IrAC complex photopleth signals. The 
magnitude of each division carries along with it the conventional (RdAC/IrAC) ratio, in 
addition, its phase encodes the angle variation between the Red and Infrared signals, which as 
disclosed in the foregoing, can generally be equal to zero for undisturbed signals. Different 



MLABS.021A 



-10- 



Knobbe, Martens, Olson & Bear, LLP 



plots in the complex plane result depending upon whether type A or type B complex 
waveforms are used to generate the complex ratios. The different plots also enable or suggest 
different types of signal processing, examples of which are disclosed as follows. 

[0056] For example, when the arterial saturation is constant and type A complex 
waveforms are generated, the complex ratio plot in the complex plane looks like a fuzzy 
point very close to the real axis, as shown in Figure 12. This is remarkable for it indicates, 
unlike previous techniques, that the instantaneous ratios over one or several photopleths are 
nearly constant throughout, which matches the input signal data, i.e., that the saturation is 
constant. Figure 19 illustrates a plot of the instantaneous saturation vs. time of the data used 
to generate the complex ratios shown in Figure 12, as well as the corresponding photopleth 
from which the saturation was calculated. Note that the saturation is generally calculated 
from the magnitudes of the complex ratios vs. time. Also, Figure 20 illustrates that the 
maximum deviation from the mean value of the saturation is less than about 0.5 percent (%), 
which compares favorably to the more than about 50 percent (%) variations calculated using 
the Biox algorithm. Moreover, Figure 20 illustrates that the standard deviation of the 
saturation is z. mere about 0.18 percent (%). These results clearly and advantageously 
indicate that weighing or filtering is not needed to utilize the data, and that all points in the 
photopleth can be useful in the subsequent analysis. 

[0057] Figure 13 is a plot of the complex ratios generated by the photopleths used 
to generate Figure 12, but with type B complex waveforms. As shown in Figure 13, the 
constant saturation translated into a circle of fixed radius in the complex plane. A myriad of 
mapping techniques available in the field of complex analysis can be brought to bear to help 
analyze this type of signal. For example, the.logarithm function may be used to map a circle 
in the complex plane into a line in the same complex plane. Figure 14 depicts a 3D plot of 
complex ratios vs. time, where time is plotted along the vertical 'Z' axis. As shown in Figure 
14, a constant ratio, i.e. constant saturation, is reflected as a uniform helix. 

Generating Confidence Measures 

[0058] Normally, each of the instantaneous complex ratio values are valid and 
reliable. However, certain physiological or patient motion conditions may render a some or 



MLABS.021A 



-11- 



Knobbe, Martens, Olson & Bear, LLP 



all of the calculated ratios unusable. Therefore, it is desirable to provide confidence 
measures by which such unreliable points may be discarded, filtered or corrected. Toward 
that end, Figure 15 illustrates a complex correlation generator 1500, which includes one or 
more conjugate generators 1502, one or more signal multipliers 1504, signal adders 1506 and 
signal dividers 1508. As shown in Figure 15, the complex RdAC and the complex IrAC are 
input into the generator 1502. The complex RdAC is multiplied by the conjugate of the 
complex IrAC. This product is divided by the sum of the complex IrAC multiplied with its 
conjugate, and the product of the complex RdAC and its conjugate. The resulting signal is 
then multiplied by the scalar two (2) to form a measure of the complex correlation, such as, 
for example, a complex confidence number associated with each calculated complex ratio. 

[0059] This confidence may be used to gate or preferentially weigh each 
corresponding complex ratio in a filtering technique to provide more reliable saturation 
values. For example, Figure 16 depicts a plot of instantaneous complex ratios calculated over 
several seconds and Figure 17 depicts the instantaneous complex correlations associated with 
the calculated complex ratios. Figure 18 depicts the joint relationship between the calculated 
complex ratios of Figure 16 and associated the complex correlation of Figure 17. Note that a 
certain correlation threshold may be established below which all data point can be rejected, 
thereby advantageously enhancing the final saturation estimation. Note that in Figure 18, all 
the data points are acceptable. Situations where not all data points are acceptable, e.g. 
abnormal waveforms or motion artifacts, will be disclosed in the following. 

[0060] The phase of the instantaneous complex ratios may also be used to assign 
a confidence measure to each complex ratio value. Typically the phase value should be very 
close to zero in the case of complex ratios generated from type A complex waveforms, as 
shown in Figure 12. However, example of the use of the phase in the filtering of unreliable 
complex ratios under the effects of motion artifacts will also be discussed in the following. 

Implementing Filtering Techniques 

[0061] Once the complex waveforms or their corresponding instantaneous 
complex ratios, or saturation values are available, a myriad of linear, nonlinear and statistical 
filtering techniques may be applied to reliably estimate the blood saturation values. For 
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example, when simple averaging is desired, the areas of the complex photopleth waveforms 
shown in Figure 6 can be calculated over a certain span of time or integral number of pulses. 
The ratio between two values corresponding to the Red and Infrared waveforms' areas is 
calculated, the result of which can then be used to calculate the blood saturation value. Note 
that this is the analog of integrating the area of the Red and Infrared photopleths shown in 
Figure 2 then taking their ratios. Figure 3 shows an example of area integration of one real 
waveform. For the case of a complex waveform the area can be estimated by summing the 
individual areas of triangles that constitute the complex waveform, as shown in Figure 7. 
Introductory calculus textbooks may be consulted on the subject of area integration. This 
technique does not require the use of complex ratios, rather straightforward real number 
division of waveform areas can be used. 

[0062] Another powerful filtering technique takes advantage of the abundance of 
ratios values available over short period of time. For example, the continuous stream of 
instantaneous ratio or saturation values can be fed into a weighing filter along with their 
associated confidence values. The filter can then discard or appropriately weigh the 
corresponding value of the ratio or the saturation. The high number of values, e.g., 62.5 
values per second in the present embodiment, available to the filter makes it more likely that 
some of them will be within an acceptable limit despite the affect of disturbances or noise. 

[0063] Statistical techniques such as frequency distribution analysis can further be 
used alone or in combination with the previous techniques to estimate the true blood 
saturation values. Exemplary techniques where the statistics of the distribution of the ratio or 
saturation values can be used to extract a more accurate estimation of the true saturation 
value, as disclosed in the following. 

Managing Abnormal Waveforms 

[0064] Under certain physiological conditions, venous blood in the vascular bed 
may undergo pulsation that may or may not be coupled to the arterial pulsation. These 
pulsation can be strong enough as to disrupt the normal ratio calculations thus giving 
erroneous saturation readings. When the arterial pulsation has a distinct frequency from the 
venous pulsation, then the arterial pulsation can be isolated using frequency analysis such as 
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the Fast Fourier Transform used to perform a Saturation Transform, as disclosed in U.S. 
Patent No 6,157,850, mentioned in the foregoing. Strong venous coupling may not be 
necessarily uniform in time across each arterial pulse, and under certain patient conditions the 
venous pulse may affect only a portion on the photopleth. This can be advantageous for an 
algorithm that analyzes the information in the time domain. Figures 21 and 22 depict 
Infrared and Red photopleths that are modulated by venous pulsation. Although their shape 
may look like a normal photopleth, the instantaneous saturation vs. time curve in Figure 23 
illustrates a large spread in the instantaneous saturation values over one photopleth, as 
compared to that of Figure 19 and Figure 20. 

[0065] Under such conditions, taking the average value (or determining the area) 
of all the instantaneous saturation points may not be the best estimate of the true saturation 
value. For exzunple, the distribution of the instantaneous saturation values of Figure 23, e.g. 
the histogram illustrated in Figure 24 shows a skewed distribution with a mode at about 
83.5 percent (%). This is about 1.8 percent (%) below the average saturation value of about 
85.3 percent (%). Another possibility is to phase filter the data and the use the distribution in 
a similar manner. Phase filtering will be discussed in more detail in the section on managing 
motion artifacts. In yet another approach, averaging the saturation values over a certain 
subsection of the photopleth, away from the location of the venous pulsation, can generate 
more accurate saturation values. 

Managing Motion Artifacts 

[0066] ; In the context of pulse oximetry, "motion artifacts 1 ' refer to any extraneous 
disturbing source that affects the shape or quality of the optically detected photopleth signal. 
The disturbance may be a deformation of the vascular bed, a decoupling of the photo detector 
from the skin, the movement of the sensor itself along the skin surface, or the like. The wide 
dynamic range of the effect, in terms of its frequency and size, as well as its multiple sources, 
makes the impact of motion artifacts on the photopleth quite dramatic. More importantly, the 
impact of motion artifacts on the calculated saturation values can also be very large, thus 
causing drawbacks in many older generations of pulse oximeters. For example, when motion 
is repetitive and affects the photopleth over its entirety, then a combination of adaptive 
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filtering and frequency-domain techniques can provide the best estimates of the saturation 
values. On the other hand, when motion is intermittent or non-repetitive, e.g. pseudo 
random, then a combination of a time-domain analysis and adaptive filtering techniques with 
fast adaptation rates works better. While each of the foregoing techniques has its strength 
and weaknesses, use of multiple parallel engines executing two or more of the foregoing 
techniques in parallel, and then fusing their results together, often provides best overall 
performance. Toward that end, the present technique of complex analysis can be a valuable 
addition to a set of parallel engines that advantageously improves the accuracy of pulse 
oximeters by correcting for a subset of conditions where the previous algorithms have failed. 
For example, the ability of the present algorithm to explicitly encode the phase of the signal 
on a point by point basis renders it very valuable in case of sensor decoupling from the skin 
where phase decorrelation between Red and Infrared photopleths is prevalent. 

[0067] Figure 25 and 26 illustrate photopleth signals that are corrupted by motion 
artifacts. Figure 27 illustrates a polar plot for type A complex waveforms generated using the 
photopleths of Figures 25 and 26, while and Figure 28 illustrates a polar plot for type B 
complex waveforms generated using the photopleths of Figures 25 and 26. It is noteworthy 
to compare and contrast Figure 12 with Figure 27. As disclosed in the foregoing, Figure 12 
illustrates complex ratios where the input signals are motion artifact free. Figure 12 exhibits 
a localized point with a magnitude of about 0.53, which generally corresponds to a saturation 
of about 99 percent (%). Note that the small angle it subtends near the real axis. On the 
other hand, Figure 27 illustrates complex ratios where the input signals are riddled with 
motion artifacts, resulting in ratios with varying magnitudes having widely varying angles. 

[0068] At the outset, it is difficult to tell which ratios of Figure 27 are true and 
which ratios were affected by noise and are therefore, false. However, as shown in Figure 28, 
a phase filter can select or pass values similar to those expected, such as, for example, values 
that subtend low angles from the origin. In one embodiment, the phase filter passes values 
corresponding to the type and value of phase angles determined, for example, through 
calibration processes performed and associated with valid data. In an embodiment, the phase 
filter selects or passes values corresponding to phase angles ranging from about -2.0 to about 
3.0 degrees, and more preferably, selects values corresponding to phase angles ranging from 
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about 0.0 to about 1.0 degrees. When the phase filtered saturation points are histogramed, as 
shown in Figure 30, the most likely saturation is somewhere near about 98 percent (%), 
which is only about 1 percent (%) away from the true saturation value of about 99 percent 
(%). Figure 31 also includes plots of unfiltered instantaneous saturation values vs. time as 
well as the phase filtered saturation values plotted on top of them. The corresponding 
photopleth is also shown for comparison. 

Using Frequency Domain Complex Ratios 

[0069] When a real signal is transformed into the frequency domain, using the 
Fourier transform for example, the corresponding frequency representation is a series of 
complex numbers. These complex numbers denote complex frequencies each having a 
magnitude and a phase. When the Red and Infrared photopleths are transformed into the 
frequency domain, their corresponding complex frequencies can advantageously be divided 
to generate complex ratios. Each complex ratio posses a magnitude and a phase similar to 
the complex ratios generated in the time domain, as disclosed in the foregoing. The 
frequency domain complex ratios are a representation that has complex ratios vs. frequency, 
as opposed to representations that have time domain complex ratios vs. time. As disclosed, 
each technique has its advantages and disadvantages depending on the type of signals present 
and the nature of the perturbations. 

[0070] An example of the use of frequency domain complex ratios according to 
aspects of the present invention are illustrated in Figure 32 and Figure 33. A series of ratios 
are calculated for the fundamental and first two (2) harmonics of the photopleths waveforms 
shown in Figure 2. Figure 32 illustrates the magnitudes of the complex frequency ratios 
calculated from the fundamental and harmonics of the photopleths of Figure 2, while 
Figure 33 illustrates frequency transformed photopleths and the phases of the corresponding 
ratios. Figure 34 illustrates the polar plot of the complex frequency ratios of Figure 32 and 
Figure 33. 

[0071] Although the foregoing disclosure generally references various signal 
processing mechanisms, a skilled artisan will recognize from the disclosure herein that the 
generators 500, 800, 1 100, and 1500 can be implemented with software, firmware, or the like 
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executing on hardware, such as, for example a digital signal processor (DSP). Moreover, the 
calculation's incorporated into the generators can be carried out using software, hardware, or 
combinations of the same. In addition, the DSP can be part of a portable or stationary device, 
such as an oximeter, personal monitoring device, or the like. 

[0072] Although the foregoing invention has been described in terms of certain 
preferred embodiments, other embodiments will be apparent to those of ordinary skill in the 
art from the disclosure herein. Additionally, other combinations, omissions, substitutions and 
modifications will be apparent to the skilled artisan in view of the disclosure herein. 
Accordingly, the present invention is not intended to be limited by the reaction of the 
preferred embodiments, but is to be defined by reference to the appended claims. 

[0073] Additionally, all publications, patents, and patent applications mentioned 
in this specification are herein incorporated by reference to the same extent as if each 
individual publication, patent, or patent application was specifically and individually 
indicated to be incorporated by reference. 
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